Energy dependence of particle multiplicities in central Au+Au collisions 
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We present the first measurement of the pseudorapidity 
density of primary charged particles in Au-|-Au collisions 
v^s^Tiv" ~ 200 GeV. For the 6% most central collisions, 
we obtain dNch/ dr]\\r,\<i = 650 ± 35(syst). Compared to 
collisions at ^s^^ — 130 GeV, the highest energy studied 
previously, an increase by a factor of 1.14 ±0.05 is found. The 
energy dependence of the pseudorapidity density is discussed 
in comparison with data from proton-induced collisions and 
theoretical predictions. 

PACS numbers: 25.75.-q 

Collisions of gold nuclei at an energy of y/s^ = 
200 GeV have been studied using the PHOBOS de- 
tector. PHOBOS is one of the experiments at the 
Relativistic Heavy-Ion Collider (RHIC) at Brookhaven 
National Laboratory aimed at understanding the behav- 
ior of strongly interacting matter at high temperature 
and density. Quantum chromodynamics (QCD), the 
fundamental theory of strong interactions, predicts that 
under these conditions, which may be probed in heavy- 
ion collisions, a new state of matter will be formed, the 
quark-gluon plasma (QGP) In this state, quarks and 
gluons are no longer confined inside hadrons, as is the 
case for normal nuclear matter. Information about the 
particle and energy density achieved in the early stages 
of the collision process is carried by the pseudorapidity 
density of particles emitted from the primary collision 
point In this analysis, we have determined the pseu- 
dorapidity density of charged particles, dNch/dr], in the 
most central Au-|-Au collisions. We focused in particular 
on the region near ij = 0, where rj = — lntan(6'/2) and 6 
is the polar angle from the beam axis. 

In combination with results from lower energies, these 
data permit a systematic analysis of particle production 
mechanisms in nucleus-nucleus collisions. Extension of 



the energy range to ^/s^ = 200 GeV allows a study of 
the relative contributions of hard parton-parton scatter- 
ing processes, which can be calculated using perturbative 
QCD, and soft processes, which are treated by phe- 
nomenological models that describe the non-perturbative 
sector of QCD. With increasing collision energy, hard 
processes are expected to become increasingly important 
for particle production near mid-rapidity. 

For Au-|-Au collisions at RHIC energies, the yield and 
momentum distribution of particles produced by hard 
scattering processes may be modified by "jet quenching" , 
i.e. the energy loss of high momentum partons in the 
nuclear medium |3j . This phenomenon has been proposed 
as a diagnostic tool for characterizing the initial parton 
density in Au-|-Au collisions at these energies. Prelimi- 
nary results indicate that in central Au-fAu collisions at 
y/s^ = 130 GeV the particle spectra at large transverse 
momenta, normalized to pp collisions, indeed change in 
comparison with lower collision energies Q . 

Early predictions for the charged particle pseudorapid- 
ity density varied by more than a factor of two, as shown 
in 1^. Data on the primary charged particle density 
dNch/dri\\n\<i at energies of ^s„„ = 56 and 130 GeV ||] 
have been analyzed in a wide variety of theoretical models 
iQ-Q. Generally, most models allowed a reasonable 
description of the energy dependence of dNch/dri\\ji\^i 
up to y^s^ — 130 GeV, with suitable choices of pa- 
rameters. For the ratio of particle densities near rj = 
at 130 GeV and 200 GeV, most calculations predict an 
increase between 9 and 20%. An interesting exception is 
the model of Wang and Gyulassy 0. When the effects 
of jet quenching are included in their calculation, they 
not only predict a suppression of particle spectra at high 
transverse momenta, but also a change of the energy 
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dependence of dNch/ d'n\\ri\<ii leading to an increase of 
more than 30% between 130 and 200 GeV when the 
default parameter set is used. The energy dependence 
of particle production presented here provides important 
constraints on this effect. 

Details of the PHOBOS experimental setup can be 
found elsewhere 16|. The apparatus, shown schemat- 
ically in Fig. 1, employs silicon detectors for vertex 
finding, particle tracking and multiplicity measurements. 
This analysis is based on the first 6 layers of the 16 
layer two-arm spectrometer (SPEC), the two-layer vertex 
detector (VTX), the single-layer octagon barrel detector 
(OCT) and the three single-layer ring detectors (RING) 
located on either side of the interaction point. The ac- 
ceptance of SPEC, VTX and OCT includes -1< < 1, 
covering different regions of azimuth. The combined ac- 
ceptance of OCT and RING detectors reaches rj = ±5.4. 
For the 2001 run period, 137168 silicon channels were 
read out, of which less than 2% were non- functional. 

The detector setup also included two sets of 16 scin- 
tillator counters ("paddle counters") located at —3.21 m 
(PN) and 3.21 m (PP) relative to the nominal interaction 
point along the beam (z) axis. These counters covered 
pseudorapidities between i < \ri\ < 4.5 and served as 
the primary event trigger. In combination with the zero- 
degree calorimeters at z — ±18.5 m, which measured the 
energy deposited by spectator neutrons, PP and PN were 
also used for offline event selection. 

Monte Carlo (MC) simulations of the detector per- 
formance were based on the HIJING event generator 
and the GEANT 3.21 simulation package, folding 
in the signal response for scintillator counters and silicon 
sensors. In the calibration of the analysis methods we 
employed a particle selection procedure that allowed us 
to iteratively modify the MC output and optimize the 
agreement of the deposited energy distributions in data 
and MC. 
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FIG. 1. Detector setup for the 2001 running period. Shown 
are the active elements of the detector systems used for this 
analysis. For clarity, one of the paddle counters and the 
two outer ring counters on each side were omitted and the 
positions of ring and paddle counters are not to scale. 



Details of the trigger performance and the event se- 
lection procedure can be found in [|6[p^. Based on MC 
simulations, we estimate that the trigger was sensitive to 
97 ± 3% of the total hadronic inelastic cross section. For 
the selection of central events used here, no contamina- 
tion from beam-gas or other background collisions was 
found. The distribution of the combined signal from PP 
and PN ("paddle mean") is shown in Fig. 2 for events 
with more than two hits in each paddle counter. This 
selection includes 88 ± 3% of the hadronic cross section. 
We selected the top 6% of the total hadronic cross section 
with the largest paddle mean, corresponding to the most 
central collisions with the largest number of participating 
nucleons. Applying the 6% centrality cut to MC events, 
we deduce that the average number of participating 
nucleons in the data is {Npart) = 344 ± 10 (syst). The 
systematic uncertainty was estimated using MC studies 
|l8t and an analysis of the correlations between the 
signals in PN and PP. 

For the selected central events, the event vertex was de- 
termined using the SPEC and VTX subdetectors. Details 
of the vertex finding procedure can be found in P,p^. 
MC studies and correlations between the vertex positions 
found independently using SPEC and VTX detectors 
indicate a vertex resolution of better than 400/im and 
a vertex finding efficiency of 100% in a fiducial range of 
— 10 cm < z^tx < 10 cm. A total of 639 central events 
were selected in this vertex range. 




2000 

Paddle mean (a.u.) 

FIG. 2. Paddle signal distribution for Au+Au events at 
^/Sj^^ = 200 GeV. The shaded area indicates the region 
selected by the 6% centrality cut. 

The charged particle pseudorapidity density was de- 
termined using four independent methods, which have 
been discussed in detail in previous publications |^,|l8|,|l9| . 
We used an analysis of two-hit combinations (track- 
lets) in SPEC and VTX detectors ||Jl|l, a hit-counting 
method in OCT and RING ||l^ and a method based on 
the observed energy loss in the silicon sensors ("analog 
method" ) ||l^ for all subdetectors. For all four methods 
the multiplicity densities dNch/dri\\jj\^i were corrected 
for particles which stop in the material of the beam pipe 
and the first detector layer, particles from secondary 
interactions and feed-down products from weak decays 
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of neutral strange particles. The resulting systematic 
errors on dA^ch/'^'7l|r;|<i for each method at ^/s'^ = 
130 GeV are described in and range from 

4.5% for the SPEC tracklet analysis to 10% for the 
hit-counting and analog methods. MC studies showed 
that the systematic uncertainties are essentially the same 
^/s^ = 200 GeV. The combined result of all four 
methods is obtained using the inverse square of the esti- 
mated uncertainties as weights in the average. Based on 
MC studies and the comparison of all four independent 
multiplicity analyses, we estimate the overall systematic 
uncertainty of the combined result to be less than 6%. 
The statistical error is negligible. We obtain a primary 
charged particle density of dA^c/i/c^?? 1 1 j,|<i = 650±35(syst) 
for the 6% most central Au-t-Au collisions at ^s„„ = 
200 GeV. Normalizing per participant pair, we find 
dNch/drj\\^l<^/{^Npart) = 3.78±0.25(syst). 
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FIG. 3. Summary of RHIC results on the pseudorapidity 
density normalized per participant pair for central Au-|-Au 
collisions at = 56, 130 and 200 GeV. The line and 

grey area show the averaged result from all 4 experiments at 



= 130 GeV. 



A compilation of all results for dNch/ drj\\j^\^i/ {^Npart) 
in central Au+Au collisions obtained at RlllC is shown 
in Fig 3. It includes results from PHOBOS at 
^s„„ — 56 GeV Q, from all four RHIC experiments 
at = 130 GeV [|jll-|2| and our new result 

Y^s„]v " 200 GeV. The average value at 130 GeV 
is dNch/dri\\,-i\^il {^Npart) — 3.37 ± 0.12, compared to 
3.78 ± 0.25 reported here for 200 GeV. 

In Fig. ^ the normalized yield per participant is com- 
pared to central Au-|-Au and Pb-|-Pb collisions at lower 
energies, including fixed target experiments [p3|-p5|. 
Also shown for comparison are results from proton- 
antiproton {pp) collisions and an interpolation of the 
pp data. The dNch/d-q values from E866/E917 ^ and 
NA49 1^,^ were obtained by integrating the measured 
charged particle rapidity and pr-distributions. Within 
the precision of the existing data, an approximately 
logarithmic rise of <iA^c/i/rf?7||,,|<i/ {\Npart) with ^s„„ is 
observed over the full range of collision energies. 



A direct measurement of the ratio R- 
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y/s^ — 200 GeV and 130 GeV was also performed. 
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FIG. 4. Energy dependence of the pseudorapidity density 
normalized per participant pair for central nucleus-nucleus 
collisions. The data are compared with pp data and nu- 
cleus-nucleus data from lower energies 
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FIG. 5. Ratio -R200/130 of pseudorapidity densities 
diVch/d77||„l<i at = 200 GeV and 130 GeV. The 

predicted increase from various models is compared to the 
observed increase in the data (vertical line). The grey band 
shows the 90% confidence interval. 

The ratio was determined separately for each of the 
four multiplicity analyses at 130 and 200 GeV and 
then averaged, weighted with the inverse square of the 
systematic error for each ratio measurement. Many 
of the systematic uncertainties in each of the analyses 
partially cancel, leading to a significant reduction in the 
uncertainty of the ratio. The theoretical uncertainty 



in the Npart calculation is also largely eliminated in 
the ratio. We estimate that the uncertainty from the 
energy dependence of absorption, feed-down and accep- 
tance corrections is less than 0.03 in i?200/i30- Based 
on MC tests and the comparison of the results for the 
four different methods we estimate a total systematic 
uncertainty of less than 0.05 for the averaged result. The 
combined estimate of the ratio of normalized multiplicity 
of nor- densities dNch/ dT]\\r,\<i/ {\Npart) at ^/s^ = 200 GeV 
and 130 GeV, obtained from all four multiplicity analysis 
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methods, is i?200/i30 = 1.14 ±0.05. The quoted uncer- 
tainty is entirely systematic and corresponds to a 90% 
confidence level. In Fig. 5 this value is compared with the 
pp parametrization |26|] and several model calculations. 
The expected increase from the interpolation oipp results 
falls slightly below the allowed region. Within the sys- 
tematic uncertainty, the result is in agreement with most 
other model predictions. However, we do not observe the 
large increase expected in the HIJING model [Q, which 
predicted i?200/i30 = 1-32 using the default jet quenching 
parameters. Further theoretical work will be required to 
describe both the preliminary results on large px spectra, 
as well as the energy dependence of particle multiplicities. 



increase in initial energy density. Over this energy 
interval, a smooth evolution of dNch/drj in all regions 
of pseudorapidity is seen. These results give further 
constraints for models including effects of the partonic 
medium in the early collision stages. 
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. The lower plot shows the shape of dNch/dr) at 
— 200 GeV, obtained from the analog method, nor- 
malized to dNch/dri\n=o- On top we show the pseudorapidity 
dependence of -R200/130 using the result of the analog method 
at both energies. In both plots, the error bars indicate the 
systematic uncertainty. 

Additional information about the change in particle 
production from 130 GeV to 200 GeV collision energy 
can be gained from Fig. 6. Using the results from the 
analog method ||T^, it shows the scaled dNch/drj out to 
ry = 5 for central collisions at 200 GeV (bottom) and 
the ratio of the pseudorapidity distributions at 200 and 
130 GeV (top). A slow increase of the ratio is observed 
to 77 « 3.5. This is followed by a steeper rise at higher 
77, corresponding to a widening of the pseudorapidity 
distribution and reflecting the shift in the fragmentation 
region due to the larger beam rapidities at 200 GeV. 

In summary, we have performed the first measurement 
of the pseudorapidity density of primary charged parti- 
cles in Au+Au collisions at y^s^ = 200 GeV. Near mid- 
rapidity, this density shows an approximately logarithmic 
evolution over a broad range of collision energies. Going 
from = 130 GeV to 200 GeV, the observed increase 

is a factor of 1.14 ± 0.05, corresponding to a moderate 
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